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IRAS 13197-1627 has them all: Compton-thin absorption, 
photo-ionized gas, thermal plasmas, and a broad Fe line 
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ABSTRACT 

We report results from the XMM-Newton observation of IRAS 13197-1627, a luminous IR 
galaxy with a Seyfert 1.8 nucleus. The hard X-ray spectrum is steep (r ~ 2.5) and is absorbed 
by Compton-thin (Nh ~ 4 x 10 23 cm' 2 ) neutral gas. We detect an Fe Ka emission line at 
6.4 keV, consistent with transmission through the absorber. The most striking result of our 
spectral analysis is the detection of a dominant X-ray reflection component and broad Fe line 
from the inner accretion disc. The reflection-dominated hard X-ray spectrum is confirmed by 
the strong Compton hump seen in a previous BeppoSAX observation and could be the sign 
that most of the primary X-rays are radiated from a compact corona (or base of the jet) within 
a few gravitational radii from the black hole. We also detect a relatively strong absorption line 
at 6.8 1 keV which, if interpreted as FeXXV resonant absorption intrinsic to the source, implies 
an outflow with velocity ~ 5 x 10 3 km s~ x . In the soft energy band, the high-resolution RGS 
and the CCD-resolution data show the presence of both photo-ionized gas and thermal plasma 
emission, the latter being most likely associated with a recent starburst of 15-2OA/ yr _1 . 
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1 INTRODUCTION 

IRAS 13197-1627 (also MCG-03-34-064, z=0.01654) is a "warm" 
(/25,im//60nm ~ 0.48) and luminous IR-galaxy with Ltr=L(8- 
1000 /xm)=1.7 x 10 n L Q (Sanders et al. 2003). It was first clas- 
sified as a moderately reddened Seyfert 2 by Osterbrock and de 
Robertis (1985) & de Robertis, Hutchings & Pitts (1998). However, 
as pointed out by these authors, the source has peculiar emission- 
line properties: in particular its Ha, [Oi], [Nn], and [Sn] are both 
exceptionally broad for Seyfert 2s (with FWHMs from 480 km s - 
to 860 km s _1 ) and asymmetric, with excess flux blue-wards of 
their centres. In the UV, most of the lines are broad at the level 
seen in the Optical, and Civ as a width comparable to those typ- 
ically measured in Seyfert 1 galaxies. Aguero et al (1994) were 
able to disentangle broad and narrow components of the Ha and 
H/9 lines which led to a Seyfert 1.8 classification. A broad compo- 
nent to Ha was also detected by Young et al (1996) in polarised 
light. Signatures of the presence of Wolf-Rayet stars have been re- 
ported by Cid Fernandes et al (2004) who also estimate that about 
25 per cent of the stellar population is relatively young (less than 
25 Myr). Extended Radio emission is also detected with a linear ex- 
tension of about 280 pc, almost aligned with the host galaxy major 
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axis (Schmitt et al. 2001a; 2001b) while the 1.4 GHz luminosity is 
1.6 x 10 30 erg s _1 Hz -1 (Condon et al. 1996). 



In the X-rays, the source was first observed with ASCA in 
1995 revealing a very steep spectral shape (photon index T ~ 3), a 
large absorbing column density TVh — 6-^8 x 10 23 cm~ 2 , a narrow 
Fe emission line at 6.4 keV, and a "soft excess" component emerg- 
ing below about 2 keV (Ueno 1997). A BeppoSAX observation was 
obtained in 1998 and was first published by Risaliti (2002) in a sta- 
tistical study of 20 bright Seyfert 2 galaxies. A variation in the ab- 
sorbing column density between the ASCA and the BeppoSAX ob- 
servations was later suggested by Risaliti, Elvis & Nicastro (2002). 
A more detailed analysis of the BeppoSAX data was presented by 
Dadina & Cappi (2004). The absorbing column was found to be 
dependent on the adopted spectral model so that the column den- 
sity variation could not be confirmed (we shall discuss this point in 
our analysis below). The better quality, broadband BeppoSAX data 
also revealed spectral complexity that could be modelled in terms 
of a partial covering scenario or of a reflection-dominated model in 
which X-ray reflection from the accretion disc dominates the hard 
spectral shape. In both cases, Dadina & Cappi estimated an intrinsic 
2-10 keV luminosity of ~ 2 x 10 44 erg s~\ making IRAS 13197- 
1627 the nearest and brightest type-1.8 quasar known to date. An 
absorption line at ~ 7.5 keV was also detected and, if interpreted as 
He-like Fe resonant absorption, the line energy implies an outflow 
with a velocity of the order of 0.1 c. 
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Figure 1. We show a soft (S: 0.3-3 keV) and ahard (H: 3-12 keV) EPIC-pn 
source light curve together with the resulting hardness ratio (bottom panel). 
As revealed by fitting a constant to the data, the hard band is more vari- 
able than the soft one (x% = 1-65 vs. Xr = 0.85 for 67 dof). However, 
the hardness ratio is still consistent with a constant value throughout the 
observation. 

Here we report on results from a new XMM-Newton observa- 
tion of IRAS 13 197-1627 and we present detailed spectral analysis 
at both CCD and gratings resolution allowing us to better under- 
stand the nature of the hard and soft X-ray emission components 
in this luminous IR galaxy. 



2 THE XMM-NEWTON OBSERVATION AND 
FIRST-LOOK ANALYSIS 

IRAS 13197-1627 was observed by XMM-Newton on 2005 Jan- 
uary 24 for a total exposure of 45 ks. The data have been processed 
starting from the ODF files with the XMM-Newton SAS software 
(version 6.5.0). Source spectra and light curves of the EPIC cam- 
eras were extracted from circular regions centred on the source, 
while background products were extracted from off-set regions 
close to the source. We remind here that the background region 
of the pn camera has to be taken as close to the centre as pos- 
sible to avoid contamination from the strong spatially dependent 
fluorescent Cu Ka emission line originating from the electronics 
board mounted below the pn CCD which would result in a spuri- 
ous absorption line around 8 keV in the (background subtracted) 
source spectrum. With our choice of the background extraction re- 
gion, such contamination is negligible. After filtering for periods 
of high background the net exposure is 37 ks in the pn camera and 
39 ks (43 ks) in the MOS 1 (MOS 2) detectors. The EPIC spectra 
were grouped to have at least 20 counts per bin to ensure the valid- 
ity of x 2 statistical analysis. The RGS were operated in the standard 
spectroscopy mode and standard data reduction was performed re- 
sulting in a net exposure of about 43 ks in the RGS detectors. 

In Fig. 1 we show the background-subtracted source light 
curve in a soft (S: 0.3-3 keV) and a hard (H: 3-12 keV) energy 
band after removal of residual background flares. No strong vari- 
ability is seen in the soft band and a fit with a constant is very satis- 
factory (xr — 0.85 for 67 dof). More fluctuations are present in the 
hard band (Xr = 1-65). In the bottom panel, we plot the hardness 
ratio and compare it with the best-fitting constant value producing 
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Figure 2. The broadband (0.3-12 keV) EPIC-pn spectrum of IRAS 13197- 
1627. The spectrum exhibits strong curvature in the hard band which is a 
clear sign of absorption by a substantial column of gas in the line of sight. 
Fe emission is detected at 6.4 keV in the rest-frame (vertical dotted line) 
together with a deep absorption edge at 7 keV. A soft excess is present below 
~ 3 keV and is characterised by a bump around 0.9 keV. Data have been 
rebinned for visual clarity. 

an acceptable fit with Xr — l-l- Given the limited evidence for 
spectral variability, in the following we consider the time-averaged 
spectrum of the source only. The broadband EPIC-pn spectrum of 
IRAS 13197-1627 is shown in Fig. 2. The MOS data are not shown 
for visual clarity and agree very well with the pn. The main fea- 
tures of the spectrum are i) strong curvature in the hard band with 
a ~ 3^1 keV low-energy cut-off and a deep spectral drop around 
7 keV where the Fe absorption edge is expected; ii) a clear narrow 
Fe emission feature at 6.4 keV (vertical dotted line in Fig. 2); iii) a 
soft excess emerging below 2-3 keV which seems structured rather 
than smooth, especially around 0.9-1 keV. In this paper we adopt 
standard cosmology parameters (Ho = 70 km s _1 , Ao = 0.73, 
and qo = 0). 



3 THE 2-12 KEV SPECTRUM OF IRAS 13197-1627 

In our analysis, we consider joint fits to the EPIC-pn and MOS (1 
and 2) cameras. The pn data are considered up to 12 keV, while 
for the MOS we limit our analysis up to 9.5 keV due to the lack of 
signal-to-noise above that energy. As a first attempt to describe the 
hard spectrum above 2 keV, we consider a simple model in which a 
power law is absorbed by neutral matter at the redshift of the source 
(the ZWABS model in XSPEC ). We also include a second power law 
component emerging at soft energies to account for the soft excess. 
The two power law slopes (Th and T s respectively) are allowed to 
be different. The observed Fe line at 6.4 keV is modelled with a 
Gaussian emission line with width, energy, and normalisation free 
to vary. All subsequent fits include Galactic absorption with column 
density fixed at its nominal value (5.8 x 10 20 cm -2 , Dickey & 
Lockman 1990). 

With this simple model, we obtain an acceptable fit to the 2- 
12 keV data with X 2 = 519 for 447 de grees of freedom (dof). The 
hard power law has a slope = 2.9 ± 0.4, marginally consistent 
with the soft one (T s = 2.4 ± 0.4). The ratio of the soft to the hard 
power law normalisations is of about 1-2 per cent only. The neutral 
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Figure 3. The 2-12 keV pn spectrum is shown together with a simple 
model (solid line) comprising a power law absorbed by a column density 
of ~ 6.2 X 10 23 cm -2 of neutral matter at the redshift of the source and 
a resolved (<r ~ 90 eV) Gaussian emission line at 6.4 keV. Only the pn 
data are shown but all spectral analysis is performed simultaneously with 
the MOS data as well. Notice several residuals: absorption structures are 
seen around 6.8 keV and 7.1 keV. Moreover the model systematically un- 
derestimates the data in the 5-6 keV band and above 10 keV. In the Figure, 
data have been rebinned for visual clarity. 



absorber affects the hard band above ~3 keV and is Compton-thin 
with a column density of 6.2 ± 0.5 x 10 23 cm -2 . 

The Fe emission line is clearly detected (A^ 2 = 160 for 2 
more free parameters) and has an energy of 6.40 ± 0.02 keV indi- 
cating an origin in neutral matter. The line has an observed equiva- 
lent width (EW) of 140 ± 20 e V when computed with respect to the 
unabsorbed continuum, to be compared with the theoretical predic- 
tion of 80 eV for a line transmitted through the observed column 
density (e.g. Matt 2002). 

We measure an Fe line width a = 901*0 eV - signifi- 
cantly broader than the EPIC cameras spectral resolution (~40 eV 
in cr). The measured Fe line width corresponds to a FWHM= 
I.21I0 33 x 10 4 km/s. Under the assumption that the line emitting 
gas is gravitationally bound and that gas motion occurs in randomly 
oriented circular orbits, the line width places the emitting gas within 
2 x 10 3 gravitational radii (r g = GM/c 2 ) from the central black 
hole (Rrolik 2001). If the line, as is likely, is transmitted through 
the Compton-thin absorber, the line width would place the absorber 
much closer to the nucleus than often thought (e.g. Galactic discs, 
Maiolino & Rieke 1995; dust lanes, Malkan, Gorijn & Tarn 1998; 
Guainazzi, Matt & Perola 2005; starburst clouds, Weaver 2001; see 
also Lamastra, Perola & Matt 2006). However, a much closer loca- 
tion of the absorbing gas, consistent with the measured line width, 
has been proposed by e.g. Elvis (2000; 2004). The case for an ab- 
sorber located close to the nucleus is particularly compelling in 
NCG 1365 (Risaliti et al. 2005) in which rapid transitions between 
Compton-thin and Compton-thick absorption have been observed. 
On the other hand, as we shall discuss below, the observed Fe line 
width could be due to an unmodelled broader component to the Fe 
line profile. 

Although our simple spectral model (hereafter the "baseline 
model") provides a reasonable description of the hard spectrum, 
several residuals are present as shown in Fig. 3. Two absorption 
features are seen around 6.8 keV and 7.1 keV, while positive resid- 
uals are left at high energy (above 10 keV) and in the 5-6 keV band 
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Figure 4. Top: A Gaussian line is added to the baseline model and its en- 
ergy and normalisation are varied. The contours represent an improvement 
of Ax 2 = -1, -2.3, -4.61, -9.21 in fitting the joint pn-MOS data. The 
Ax 2 = — 1 contour (outermost) is shown as a reference of the best-fitting 
continuum model. Bottom: In this case the Fe Ka line of the baseline model 
is forced to be unresolved. The shape of the residuals strongly suggests 
the presence of X-ray reflection from the disc (broad Fe emission line and 
broad Fe absorption edge). 



where the model seems to systematically underestimate the data. 
As a first check for the significance of the residuals, we add an un- 
resolved Gaussian line to the baseline model, vary its energy in the 
4-8 keV range, its normalisation in the ±1 x 10~ 5 ph cm -1 s _1 
range, and record the \ improvement (see Miniutti & Fabian 2006 
for more details). The results are presented in the top panel of Fig. 4 
where we show the 68, 90, and 99 per cent confidence contour lev- 
els for the additional Gaussian line rest-frame energy and intensity 
in the most relevant band around the Fe K complex. The outermost 
line shows the A\ 2 = — 1 contour as reference. 

Although the contours should not be taken as a precise signif- 
icance test, they indicate the relative importance of the residuals. 
Fig. 4 clearly reveals both positive and negative residuals on the 
red and blue side respectively of the Fe Ka line (6.4 keV). If the 
Fe Ka emission line in the best-fitting model is forced to be nar- 
row (i.e. unresolved, as in most Compton-thin sources), the bottom 
panel figure is obtained, where the positive residuals are very rem- 
iniscent of X-ray reflection from the accretion disc resulting in a 
skewed and broad Fe emission line and Fe absorption edge. With 
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Table 1. The best-fitting parameters from joint fits to the 2-12 keV (pn) 
and 2-9.5 keV (MOS) spectrum of IRAS 13197-1627 (x 2 = 442 for 438 
dof). 
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Figure 5. Same as in Fig. 4, but now for the best-fitting model comprising 
Compton-thin absorption and reflection from the accretion disc, both with 
slightly super-solar Fe abundance (see text for details). The only remaining 
residual is a narrow absorption line at 6.8 keV. The rest-frame energies of 
the closest transitions (FeXXV and FeXXVI resonant absorption) are shown 
as vertical lines. 



the caveats mentioned above, the significance of the whole struc- 
ture seem to be higher than 99 per cent. A more detailed spectral 
modelling of these absorption and emission features is presented 
below. 



3.1 X-ray reflection from the accretion disc and Fexxv 
resonant absorption 

The shape of the residuals in Fig. 4 strongly suggests that both the 
positive and negative residuals could be accounted for by a reflec- 
tion component from the accretion disc. We then add a relativisti- 
cally blurred reflection spectrum in which reflection continuum and 
emission lines are computed self-consistently (from Ross & Fabian 
2005). The model free parameters are the ionization parameter f, 
the illuminating power law slope, the normalisation, and the Fe 
abundance. To test whether the Fe abundance is different from so- 
lar, we let it free to vary, replace the Compton-thin absorber model 
with one which allows to vary the Fe abundance (ZVFEABS) and 
forced the Fe abundance in the two models to be the same. The 
illuminating power law slope is forced to have the same photon 
index as the hard power law continuum. The relativistic blurring 
is obtained by using a LAOR kernel in which the emissivity pro- 
file is fixed at its standard value (e = r~ 3 ), the outer disc radius 
at 400 r g , while the inner disc radius and observer inclination are 
free parameters. The reflection spectrum is absorbed by the same 
column density as the continuum. 

We obtain a very significant improvement with respect to 
the previous model, with A\ 2 = 63 for 5 more free parameters 
(X 2 = 456 for 442 dof). Since the reflection model does not in- 
clude emission from Ni (see Ross & Fabian 2005), we add a Gaus- 
sian emission line with energy between 7.4 keV and 7.6 keV and 
apply to it the same relativistic blurring as for the reflection spec- 
trum. We obtain a marginal improvement of A\ 2 = 6 for 2 more 
free parameters. Before discussing the best-fit parameters, we con- 
sider again the procedure used to produce Fig. 4 on this best-fitting 
model to search for any other possible residual (see Fig. 5). The X- 
ray reflection from the disc accounted for all positive and negative 
residuals (compare Fig. 5 with Fig. 4) except what appears to be 



Parameter 


Value 


F-test 


Continuum 


r h 

iv^ pn cm s ] 

r s 

N s [1CT 5 phcm" 2 s" 1 ] 


2.5 ±0.3 

c 7 +1.6 
b -'-1.9 

2.1 ±0.4 

8.2 ±4.8 




Absorption 




> 99.99% 


N H [10 23 cm" 2 ] 


3.9 ±0.4 




Fe abundance 




99.8% 


Ap e [solar] 


1.5 ±0.3 




Transmitted Fe Ka line 




> 99.99% 


E [keV] 
a [eV] 
EW [eV] 


6.40 ± 0.02 
< 85 
105 ± 18 




Disc reflection 




99.8% 


£ [erg cm s~ x ] 
Em [keV] 
EW Nl [eV] 
R 

r m [r g = GM/c 2 ] 
i [degrees] 


< 60 
7.5 ±0.06 
60 ±40 
10 ±3 

27 ± 17 


94.8% 
99.1% 


Fexxv absorption line 




98.1% 


E [keV] 
EW [eV] 


fi 01 +0.08 

-50 ± 30 




Flux and Luminosity 


F2-10 [1CT 12 erg cm" 2 s" 1 ] 
L°^ 10 [10 43 ergs" 1 ] 
Lfl 10 [lO^ergs" 1 ] 


2.1 ±0.2 
0.7 ±0.2 
4 2+ 2 ' 8 





The EW of the Fe Ka emission line is computed with respect to the unab- 
sorbed continuum assuming it is transmitted and not absorbed. If the line 
is assumed to be absorbed by the same column as the continuum its EW is 
160 ± 70 eV. The F-test results are obtained by removing the relevant com- 
ponent from the model and by re-fitting the data. The F-test corresponding 
to r; n in the Disc reflection section is relative to the relativistic blurring 
kernel, i.e. it provides the significance of the relativistic effects. 

an absorption line at ~ 6.8 keV whose significance can be roughly 
estimated to be between 90 and 99 per cent. 

We then add a further narrow Gaussian absorption line to our 
spectral model and find a marginal improvement of Ax 2 = 8 for 2 
more free parameters (98.1 per cent significance level according to 
the F-test) for a line energy of 6.81^Q Qg keV with equivalent width 
of -50 ± 30 eV and a final statistics of x 2 = 442 for 438 dof. The 
absorption line energy suggests an origin in resonant absorption 
due to FeXXV or FeXXVI (see e.g. Bianchi et al. 2005b) but is in- 
consistent with the expected rest-frame energies (weighted mean of 
6.697 keV and 6.966 keV respectively). If FeXXV is assumed, the 
line energy indicates an outflow velocity of 5061^27oo km s_1 - It 
should be stressed that the inferred outflow velocity is remarkably 
similar to the systemic velocity of the galaxy. IRAS 13197-1627 
is at redshift z=0.0 16541 which corresponds to a receding velocity 
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Figure 6. In the top panel we show the result of our best-fitting model in 
the 2-12 keV band (see Table 1). Only the pn data are shown for clarity. In 
the bottom panel, we show the model components. 



of 4959 km s~ . Such a coincidence raises the possibility that the 
line has a local origin (i.e. Galactic) as it was recently suggested 
by McKernan, Yaqoob & Reynolds (2004; 2005) for the cases of 
extreme relativistic outflows inferred from X-ray spectroscopy in 
several AGNs (see e.g. Pounds et al 2003; Reeves, O'Brien, & Ward 
2003). 

The data, residuals, and best-fitting model are shown in the 
two panels of Fig. 6. We include the Ni emission line from the disc 
and the 6.8 keV absorption line in the final model, but we stress 
again that they are respectively detected at about the 95 and 98 
per cent level only. The best-fitting parameters are reported in Ta- 
ble 1 where we also compute the F-test significance for the most 
important spectral components obtained by re-fitting the data after 
having excluded the relevant component. In the case of the disc re- 
flection component, the F-test is computed for the overall blurred 
reflection (99.8%), for the addition of the Ni emission line (94.8%), 
and for the relativistic blurring only (99.1%) separately. 

The most important result is that we infer a very strong (and 
nearly neutral) reflection component from the accretion disc which 
dominates the hard spectrum. Notice also that the transmitted Fe 
line width is now only an upper limit, consistent with an absorber 
far away from the nucleus. The transmitted Fe line equivalent width 
is also slightly reduced because of the broad Fe line associated with 
the disc reflection component. Part of the hard curvature is now 



due to the disc reflection component which results in a slightly 
smaller column density of the absorber. As for the disc parame- 
ters, the inner disc radius is consistent with both a spinning and a 
non-spinning black hole, while the inclination is relatively poorly 
constrained. 

The large value of the reflection fraction R — 10 ± 3 (where 
R = Q/2tv and is the solid angle subtended by the reflector) is 
unusual and the resulting reflection-dominated spectrum has con- 
sequences on luminosity estimates of the nuclear emission. If the 
spectrum is indeed reflection-dominated, the intrinsic AGN lumi- 
nosity is much higher than the observed one. For this reason, in 
Table 1 we report two different measures of the 2-10 keV luminos- 
ity, with two different meanings. The first one is the absorption- 
corrected observed luminosity (Z/2- s io)> while the second is the 
estimated intrinsic luminosity of the AGN (Lfl 10 ): this is the 
absorption-corrected luminosity of the direct AGN emission (the 
hard power law continuum) multiplied by the reflection fraction 
i.e. the observed AGN luminosity plus the luminosity needed to 
produce the reflection spectrum we observe. It is clearly a model- 
dependent, but in our opinion meaningful, estimate. 

3.2 A partial covering alternative 

In their analysis of the previous BeppoSAX observation, Dadina & 
Cappi (2004) pointed out that, besides the Compton-thin absorber, 
the BeppoSAX hard spectrum could be described either in terms of 
a strong and dominant X-ray reflection component from the disc or 
in terms of a partial covering model. Their partial-covering solu- 
tion was an attempt to describe the 5-6 keV and >10 keV positive 
residuals seen in Fig 4 (and visible in the BeppoSAX data as well) 
with the curvature produced by a high-column (~ 5 x 10 24 cm -2 ) 
absorber. Motivated by their study, we also considered a similar 
model and allowed the Fe abundance of the Compton-thin absorber 
and of the partial coverer to vary, but forcing it of course to be 
the same. We also include the 6.8 keV absorption line. However, 
the partial covering best-fit is statistically worse than our best-fit 
model with x 2 = 470 for 442 dof (to be compared with \ 2 = 442 
for 438 dof). The partial-coverer has Nh — 5 x 10 24 cm -2 and 
a covering fraction of about 97 per cent, while the Compton-thin 
absorber parameters are consistent with those reported in Table 1 . 
The partial coverer accounts for the residuals above 10 keV (see 
Fig. 3) but not for the positive residuals in the 5-6 keV band and 
the absorption structure around 7 keV is only poorly modelled (see 
Fig. 4). Therefore, we do not consider the partial covering model 
any further. 



4 THE PREVIOUS ASCA AND BEPPOSAX 
OBSERVATIONS 

IRAS 13197-1627 has been observed twice previously in X-rays, 
in July 1995 by ASCA in the 0.5-10 keV band for 37 ks (Ueno 
1997), and in July 1998 by BeppoSAX in the 0.5-100 keV band 
for 44 ks, reduced to 20 ks above 10 keV (Risaliti 2002; Dadina 
& Cappi 2004). The ASCA data were retrieved from the Tartarus 
database and the event files were processed as standard and veri- 
fied to give consistent results with the products obtained by Tar- 
tarus team. We used data from the SIS and 1 and from the GIS 2 
and 3 detectors and performed joint fit to the X-ray spectra in the 
2-10 keV band. The BeppoSAX data were retrieved from the Bep- 
poSAX Science Data Centre and event files were processed as stan- 
dard. The final MECS and PDS spectra were found to be in excel- 
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Figure 7. In the top panel we show the baseline XMM-Newton model ap- 
plied to the BeppoSAX data (MECS and PDS detectors). The fit is accept- 
able but large residuals are seen around 20-30 keV where reflection is ex- 
pected to dominate. In the bottom panel, we add the disc reflection compo- 
nent, i.e. we apply the XMM-Newton best-fitting model to the BeppoSAX 
data. 

lent agreement with those from the automated pipeline (available 
on-line). 

4.1 Confirming the reflection-dominated spectrum with 

BeppoSAX 

Besides the broad Fe line and edge, the X-ray reflection model pre- 
dicts the presence of a strong Compton hump around 20-30 keV. 
The previous BeppoSAX observation can thus be used to confirm 
the detection of the Compton hump in the high-energy PDS de- 
tector. We first apply to the 2-80 keV BeppoSAX data our baseline 
model (no disc reflection). As already reported by Dadina & Cappi 
(2004) an absorption line is tentatively detected at 7.5 ± 0.2 keV, 
while none is seen at 6.8 keV. If related to resonant absorption from 
FeXXV (FeXXVl), the energy of the line implies an outflow of at 
least 27000 (13000) km s _1 . We confirm their result and add that 
the line is detected only marginally (at the 98 per cent level) but if 
true, its presence could indicate that an outflow with variable ve- 
locity on long timescales is indeed present in IRAS 13197-1627. 

The BeppoSAX data and baseline (plus absorption line) best- 
fitting model are shown in the top panel of Fig. 7. The fit is ac- 



Table 2. The best-fitting model to the XMM-Newton data is applied to 
the previous ASCA (2-10 keV) and BeppoSAX (2-80 keV) observations. 
We report the most relevant best-fitting parameters and the resulting statis- 
tics (x 2 /dof). For the BeppoSAX (XMM-Newton) observation, the baseline 
model also comprises a ~7.5 keV (~6.8 keV) absorption line. 



Parameter. ASCA 

(1995) 


BeppoSAX 

(1998) 


XMM-Newton 

(2005) 




2 6 +0 - 5 


n 9+0.7 
z - -0.4 


2.5 ±0.3 


N H 


4 3+1.3 


3.7 ± 0.8 


3.9 ±0.4 


EFe 


6.5 ±0.1 


6.4 ± 0.1 


6.40 ± 0.02 


N Fe 


6.5 ±0.9 


6.8 ±0.8 


6.2 ±0.4 


-Fcont 


9.2 ± 1.6 


11.4 ±1.3 






7.0 ±2.7 


8.0 ± 1.1 


3.9 ±0.8 


Best Fit 
Baseline 


156/150 
165/153 


52/66 
76/69 


442/438 
505/445 



The column density of the absorber (Nh) is in units of 10 23 cm -2 and 
the Fe line energy (Ep e ) is in keV. The Fe line normalisation (Np e ) is in 
units of 10~ 6 ph cm' 2 s" 1 . The continuum and reflection fluxes (F cont 
and -F rc fl) are unabsorbed and given in the 2-10 keV band in units of 
10~ 12 erg cm -2 s _1 . 

ceptable (\ 2 = 76 for 69 dof), but clear residuals are left around 
20-30 keV, where reflection is expected to dominate. We thus con- 
sider our XMM-Newton best-fitting model comprising the X-ray 
reflection component from the disc. Given the poorer quality of the 
data, we fix the Fe abundance, inner disc radius, inclination and Ni 
emission line energy to the best-fitting XMM-Newton values and 
fit the BeppoSAX data with the other parameters free to vary. The 
result is shown in the bottom panel of Fig. 7 and shows that the 
X-ray reflection component is required by the PDS data. We ob- 
tain a final fit of \ 2 = 52 for 66 dof, significantly better than the 
baseline model one. BeppoSAX measures a reflection fraction of 
R = 11 ± 2, consistent with the XMM-Newton one (R = 10 ± 3). 

4.2 Long-term spectral variability 

We have studied the long-term spectral variability of IRAS 13197— 
1627 by using the three available X-ray observations with the goal 
of determining i) changes in the continuum slope, column density 
of the absorber, and narrow Fe line flux, and ii) variability of the 
two main spectral components, namely the continuum power law 
and the disc reflection. We mention here that, based on literature 
results from ASCA (Ueno 1997) and on the Risaliti (2002) analy- 
sis of the BeppoSAX data, Risaliti, Elvis & Nicastro (2002) have 
pointed out that a long-term variability in the absorbing column 
density could be claimed in IRAS 13197-1627, which would be 
consistent with a clumpy absorbing medium close to the central 
engine. 

We applied our best-fit model to the ASCA spectrum as well 
and our results are reported in Table 2 for the three observations. 
The source observed (i.e. absorbed) 2-10 keV flux was low during 
the XMM-Newton observation (2.1 ± 0.2 x 10~ 12 erg cm -2 s" 1 ), 
about a factor 2 higher during the BeppoSAX one (4.7 ± 0.7 x 
10 -12 erg cm -2 s _1 ), and intermediate during the ASCA pointing 
(3.5 ± 0.9 x 10~ 12 erg cm" 2 s" 1 ). Table 2 shows that the hard 
power law photon index, the column density of the Compton-thin 
absorber, and the transmitted Fe line normalisation are all consis- 
tent with being the same in the three observations. We conclude that 
the data collected so far are consistent with the flux variability oc- 
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Figure 8. The 2-10 keV reflection flux is plotted against the 2-10 keV con- 
tinuum flux from the three available observations. This demonstrates that 
the reflection component responds to the continuum long-term variability 
within the errors. We also show as a reference the best-fitting linear rela- 
tionship. 

curring at constant spectral slope and absorption, meaning that we 
have no observational evidence to support the conclusion that the 
absorber is located relatively close to the central engine (although 
we cannot exclude it). 

As for the variability of the continuum and reflection compo- 
nents, we find that they are consistent with being correlated. This 
is shown in Fig. 8 where the 2-10 keV flux of the reflection com- 
ponent is plotted against the 2-10 keV flux of the continuum (both 
unabsorbed). The origin of the reflection-dominated spectrum and 
of its variability are discussed below 



5 ORIGIN OF THE REFLECTION-DOMINATED 
SPECTRUM 

The most remarkable result of our analysis is the presence of a dom- 
inant reflection component from the accretion disc with a reflection 
fraction of R = 10 ± 3. In other words, the hard X-ray spectrum 
of IRAS 13197-1627 is largely reflection-dominated. Since reflec- 
tion comes from the inner accretion disc, and since most of the 
accretion power has to be radiated from there (or in a corona above 
it), it is very difficult to imagine a situation in which the spectrum 
is reflection-dominated because absorption covers the direct con- 
tinuum X-ray source and not the inner disc. We tried anyway a 
spectral model in which the hard power law and the reflection com- 
ponent are absorbed by different column densities to check whether 
this could help explaining the large reflection fraction. However, we 
did not find any statistical improvement and, more importantly, the 
two column densities turned out to be consistent with each other 
and with the value reported in Table 1. 

One possibility is that the large observed reflection fraction is 
due to strong gravity effects. If the primary emitting source of X- 
rays is located only a few r fl from the black hole, light bending fo- 
cus the continuum towards the accretion disc (and black hole) dra- 
matically reducing the primary continuum at infinity and enhanc- 
ing the reflection fraction up to the very large observed values (see 
the light bending model proposed by Miniutti et al. 2003; Mini- 
utti & Fabian 2004; see also Fabian et al 2005; Ponti et al 2006; 
Miniutti et al 2006 for recent applications). As for the variability, 



the model was devised to reproduce large continuum variation with 
no or little reflection variability (as observed e.g. in MCG-6-30- 
15, see Fabian & Vaughan 2003; Miniutti et al 2006). The model 
successfully reproduces this behaviour for sources with R between 
~ 1 and ~ 3 — 4(suchas MCG-6-30-15), while it predicts a corre- 
lated variability for reflection-dominated sources with R > 3 — 4. 
Thus, in the present case, the continuum/reflection correlated vari- 
ability (see Fig. 8) is consistent with the model (since 7? ~ 10). 
However, a correlation between disc reflection and continuum is 
expected in any simple disc reflection model and is not specific of 
the light bending one. The main advantage of the model is that the 
large reflection fraction is naturally explained. 

As recently pointed out by Merloni et al. (2006), reflection- 
dominated spectra could also be produced from discs which are 
subject to instabilities making both density and heating rate in the 
inner disc inhomogeneous. If the clouds generated by the inho- 
mogeneous flow have high effective optical depth, steep spectra 
(r ~ 2.4) emerge and they are associated with a dominant reflec- 
tion component, in good agreement with our analysis. It is a natural 
consequence of this model that some of the inner-disc clouds re- 
sponsible for the X-ray reflection would be seen in absorption with 
a range of column densities (say 10 23 — 10 25 cm -2 ). Given the lo- 
cation of the clouds (few innermost r g ) short-timescale variability 
in absorption has to be expected and could be used to disentan- 
gle between the two models in future longer observations. At the 
present time, we have no evidence for absorption variability (no- 
tice that the Compton-thin absorber we detect has nothing to do 
with this scenario since the transmitted Fe line is too narrow to be 
produced in the inner-disc). 



6 THE RICH SOFT X-RAY SPECTRUM 

Since the source primary emission is relatively heavily obscured 
by a ~ 4 x 10 23 cm -2 column density, the AGN emission does 
not contribute directly in the soft band. This is typical of Seyfert 2 
galaxies in which the soft emission is mainly due to a blend of emis- 
sion lines for He-like and H-like K transitions of light elements 
and L transitions of Fe (e.g. Kinkhabwala et al. 2002; Bianchi et 
al. 2005a). The emission lines are most likely from gas photo- 
ionized by the AGN and located on large scales (possibly coin- 
cident with the Narrow-Line region, see e.g. Bianchi, Guainazzi 
& Chiaberge 2006) and/or collisionally ionized plasmas (e.g. asso- 
ciated with star-forming regions or gas shocked by a jet). In both 
cases a residual weak continuum due to electron-scattering of the 
primary AGN radiation by extended gas or to part of the nuclear 
emission leaking through the absorber can be present. We first con- 
centrate on the high-resolution RGS spectrum and search for signa- 
tures of emission lines related to photo-ionized gas and/or thermal 
plasmas emission. 

6.1 The high-resolution RGS data 

6.1.1 Photo-Ionized plasma 

The first order 0.45-0.7 keV RGS spectrum is shown in the top 
panel of Fig. 9 where clear Ka emission lines from N VII, O VII, 
and O VIII are seen (the O VII line is actually a He-like multiplet 
comprising forbidden, intercombination, and resonance lines). In 
all subsequent fits, a phenomenological power law continuum ab- 
sorbed by the Galactic column density is considered and its prop- 
erties will be discussed in detail in the subsequent sections. 
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required by the data at the 96 per cent level (F-test), but its slope 
and normalisation are largely unconstrained because most of the 
RGS photons are in the emission lines rather than in the contin- 
uum. The issue will be re-examined when the soft EPIC spectrum 
is discussed below. 



6.1.2 Collisionally— ionized plasma emission? 

When higher energy data up to 1.2 keV are considered, large resid- 
uals are seen and several additional emission lines are required 
above 0.7 keV (see Fig. 9, bottom). Re-fitting the XSTAR model 
discussed above in the 0.45-1.2 keV band produces an unaccept- 
able fit with % 2 = 263 for 127 dof. We then add further Gaussian 
emission lines as required in order to account for the residuals seen 
above 0.7 keV. The lines we detect are reported in Table 2 together 
with their most likely identification and with the properties of the 
power law continuum (see caption). We clearly detect Ne IX and 
Ne X Kq lines and a plethora of other emission features from Fe 
L. The Fe L emission complex (mostly from Fe XVI and Fe XVII) 
is generally interpreted as a sign of collisionally-ionized plasma 
emission. However, as pointed out e.g. by Kinkhabwala et al. 2002 
(and many others) ambiguities between photo-ionized gas and ther- 
mal plasma emission are difficult to disentangle. One possibility is 
represented by the detection of radiative recombination continua 
(RRC) or higher order series transitions for O VII (/3, 7, 8) which 
would rule out a very strong contribution from a thermal plasma 
(e.g. Kinkhabwala et al. 2002). In the present case, however, the 
statistics in the RGS is limited and we cannot easily disentangle be- 
tween the two possibilities. In the following we present a possible 
solution in terms of photo-ionized plus collisionally-ionized emis- 
sion, but firm conclusions must await better quality high-resolution 
spectra. 



Figure 9. Top: a portion of the first order RGS spectrum in the 0.45- 
0.7 keV band. The vertical lines are the theoretical energies of the emission 
lines shifted to account for the source redshift, while the solid line is the 
best-fitting model comprising a power law absorbed by the Galactic column 
plus four Gaussian emission lines. The dominant line in the He-like O VII 
multiplet is the forbidden, suggesting emission from photo-ionized plasma 
as typically observed in Seyfert 2 galaxies. Bottom: the RGS data up to 
1.2 keV are shown together with the best-fitting XSTAR model. Photo- 
ionized gas emission alone does not seem to be able to reproduce the emis- 
sion lines seen above ~0.7 keV. 



We detect Ka emission lines from N VII at 0.50 ± 0.05 keV 
and O VIII at 0.653 ± 0.006 keV, while the O VIII k/3 line is only 
an upper limit. As for the O VII He-like multiplet, we clearly de- 
tect a strong forbidden (f) line at 0.561 ± 0.002 keV and a fainter 
resonance (r) one at 0.573 ± 0.004 keV. All observed energies 
are consistent with the rest-frame theoretical energies and are re- 
ported together with the corresponding fluxes in Table 2. The O VII 
forbidden/resonance flux ratio is 2.4^g'g. The dominance of the 
O VII forbidden line relative to the resonance is generally taken as 
a signature of emission predominantly from photo-ionized rather 
than collisional plasmas (e.g. Porquet & Dubau 2000 and refer- 
ences therein) and we conclude that such a component is indeed 
present in the soft X-ray spectrum of IRAS 13197-1627. We also 
modelled the data with the photo-ionized emission code XSTAR 
and find a good fit (x 2 = 53 for 46 dof) with a column density of 
~ 1 x 10 22 cm -2 and ionization log £ ~ 1 erg cm s _1 (both loosely 
constrained). The power law continuum is (although marginally) 



6.1.3 A global model to the RGS data 

A good quality global fit to the 0.45-1.2 keV RGS data can be 
obtained by considering a composite model comprising the X-ray 
continuum (power law absorbed by the Galactic column), a photo- 
ionized gas emission component as obtained with the XSTAR code, 
and a thermal plasma component modelled through the VMEKAL 
model to account for the Fe L complex. Given that the 2-12 keV 
spectrum provides indication that Fe is overabundant, we let the Fe 
abundance in the VMEKAL and XSTAR models free to vary and 
tied to each other, while all the other elements have fixed solar 
abundances. We obtain a good global fit (x 2 = 153 for 125 de- 
grees of freedom). The higher energy lines (see Fig. 9) are well 
accounted for by the thermal model and we measure a tempera- 
ture of 0.63t ;g 4 keV, while the Fe abundance is 1.3 ± 0.3 times 
solar. The VMEKAL model improves the fit by A^ 2 = 110 for 

3 dof with respect to the power law plus XSTAR spectral model. 
We point out that we were not able to obtain a similar improve- 
ment by considering two photo-ionized components with differ- 
ent ionization parameters and column densities suggesting that a 
non photo-ionized component is indeed required by the data. Some 
positive residuals are still left around 1 keV and they can be mod- 
elled by adding a second higher-temperature thermal component. 
Its temperature is only poorly constrained in the range of 1.5- 

4 keV but the additional model improves further the statistics with 
Ax 2 = 20 for 2 more free parameters. According to our best- 
fitting model extrapolated to 2 keV, the photo-ionized component 
has a 0.5-2 keV luminosity of 5 ± 2 x 10 40 erg s~\ while the two 
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Table 3. List of the emission lines included in the best fit model for the 
RGS data between 0.45 keV and 1 .2 keV. The 90 per cent errors are given 
in parenthesis on the last digit. When required, the upper level (UL) and 
lower level (LL) of the most likely transition are reported. The model also 
comprises a poorly constrained continuum, modelled with a power law with 
T a = 2.1 ±0.8 and normalisation of 3. 5 ±3.0 X 10~ 5 phcm" 2 s _1 (we 
assume Galactic absorption). 



Eobs(keV) 


Flux" 


Likely Id. 
UL 


Etheo (keV) 
LL 


0.50(5) 


1 8(6) 


N VII Ka 


0.500 


0.561(2) 


3.8(9) 


O VII Ka 


0.561 (f) 


0.573(4) 


1.6(5) 


O VII Ka 


0.569 (i); 0.574 (r) 


0.653(6) 


2.0(5) 


O VIII Ka 


0.654 


0.666* 


< 1 


OVIIK/3 


0.666 


0.73(2) 


1.0(5) 


FeXVI 

2 P 2 l/2 2 Pl/2 3s 


0.727 

2 P?/2 2 P3/2 


0.74(2) 


1.0(7) 


FeXVI 
2 P 2 1/2 2pl /2 3s 


0.739 

2 Pl/2 2 P3/2 


0.78(2) 


0.8(4) 


Fe XVII 
2 P 2 1/2 2pj /2 3s 


0.775 

2 Pl/2 2 P3/2 


0.827(3) 


1.3(5) 


FeXVI 

2pi/22pj /2 3d 3/2 


0.826 

2 Pl/2 2 P3/2 


0.85(2) 


1.0(4) 


FeXVII 

2 P?/2 2 P3/2 M 5/2 


0.851 

2 Pl/2 2 P3/2 


0.870(3) 


1.1(5) 


Fe XVII 

2pi/22p^ /2 3fi 3/2 


0.867 

2 Pl/2 2 P3/2 


0.90(1) 


0.7(3) 


Ne IX Ka 


0.905 (f) 


0.92(1) 


1.3(5) 


Ne IX Ka 


0.915 (i); 0.922 (r) 


1.02(2) 


0.9(6) 


NeX Ka 


1.022 


1.11(2) 


0.8(5) 


FeXVI 

2pi/22p3 /2 5d 3/ 2 


1.112 

2 Pl/2 2 P3/2 




Observed Energy (keV) 




0.5 1 
Observed Energy (keV) 

Figure 10. Top: the 0.45-1.2 keV EPIC pn and MOS spectra fitted with 
the RGS model comprising a power law continuum, emission from photo- 
ionized gas modelled through the XS TAR code and two thermal components 
(VMEKAL). Bottom: the best-fitting model is shown in the 0.2-2 keV band 
(grey). Besides the soft power law, the photo-ionized gas and the two ther- 
mal plasma spectral models are shown. The overall spectrum is absorbed by 
the Galactic column density. 



a Fluxes in units of 10 5 ph cm 2 s 1 . 
* Indicates that the parameter has been fixed. 

thermal components dominate the soft band with a luminosity of 
1.4 ± 0.4 x 10 41 erg s" 1 . The soft power law contribution is much 
more uncertain (8±7x 10 40 erg s _1 ) and will be better constrained 
with the higher statistics EPIC spectrum as discussed below. 

6.2 The CCD-resolution EPIC data 

We have applied our best-fit RGS model to the pn and MOS data 
in the same 0.45-1.2 keV energy band. This allows us to i) check 
the consistency between the high-resolution and CCD-resolution 
spectra, and ii) explore whether the soft power law continuum is 
truly required by the data. We obtain a good description of the 



EPIC data with a final result of \ 2 = 273 for 233 dof with a two- 
temperature model with 0.6 ± 0.2 keV and 1.4 ± 0.2 keV plus 
photo-ionized gas emission. The Fe abundance is measured from 
the XSTAR and VMEKAL models as 1.4+o'J times solar which is 
fully consistent with the Fe overabundance we found modelling 
the 2-12 keV hard spectrum (see Table 1). However, if the Fe 
abundance is forced to be solar, the worsening of the statistics is 
marginal and a super-solar abundance is required at the 95 per cent 
level only. Finally, the power law continuum is very significantly 
required by the data: if it is excluded from the model, the statistics 
is worse by A\ 2 = 32 (2 dof). We measure a slope F s = 2.3 ± 0.3 
and a normalisation of 7.0 ± 2.0 x 10~ 5 ph cm~ 2 s _1 . Our results 
are summarised in Table 4 and remain unchanged in the broad- 
band analysis discussed below. The EPIC spectra, best-fit model 
and data to model ratio are shown in the top panel of Fig. 10, while 
the best-fitting model components are shown in the bottom panel. 
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Table 4. Parameters of the best-fit model to the 0.45-1.2 keV EPIC spec- 
trum. The model comprises Galactic absorption, a power law continuum, 
two thermal components (VMEKAL) and a photo-ionized gas emission one 
(XSTAR). The Fe abundance is forced to be the same in the two thermal 
and the photo-ionization models. For each component we also report the 
luminosity in the 0.5-2 keV band (we checked these luminosities also in 
the broadband analysis discussed below and find excellent agreement). All 
components are statistically very significantly required. However, forcing a 
solar Fe abundance produce a worsening of the fit at the ~ 95% level only. 



Parameter 


Value 


F-test 


Continuum 




> 99.997o 


r s 

L s [10 41 ergs" 1 ] 


2.3 ± 0.3 
1.0 ±0.3 




VMEKAL 1 




> 99.99% 


KT\ [keV] 

L\ h [10 41 ergs" 1 ] 


0.6 ±0.2 
0.8 ±0.2 




VMEKAL 2 




> 99.99% 


KT 2 [keV] 

Lf [10 41 ergs- 1 ] 


1.4 ±0.2 
0.6 ±0.2 




XSTAR 




> 99.99% 


N H [10 22 cm" 2 ] 
log § [erg cm s — 1 ] 
L ph [1Q 4i ergs- 1 ] 


0.9 ±0.5 
0.8 ±0.4 
0.4 ±0.2 




Fe abundance 




> 94.8% 


Ap c [solar] 


1 4+0-5 
4 - 4 -0.3 





7 THE BROADBAND 0.45-12 KEV SPECTRUM 

After having analysed the X-ray spectrum of IRAS 131 97- 1 627 by 
looking separately at the different energy bands and available en- 
ergy resolutions, the task of describing the broadband EPIC spec- 
trum of the source is greatly simplified. We consider here the best- 
fit reflection-dominated model to the hard spectrum (Table 1) and 
add to it the soft model we just discussed above. The Fe abundance 
in the Compton-thin absorber, the reflection component, the two 
thermal models, and the photo-ionized gas one is forced to be the 
same for consistency. We find an excellent description of the broad- 
band 0.45-12 keV EPIC data with x 2 = 924 for 820 dof. The spec- 
tra, best-fit model, and residuals are shown in Fig. 11. The model 
components are indistinguishable from those shown in the bottom 
panels of Fig 6 (hard band) and Fig. 10 (soft band). 

The best-fitting parameters for the reflection component, rel- 
ativistic blurring, 6.8 keV absorption line, photo-ionized and col- 
lisional plasma emission are consistent with Table 1 (hard band) 
and with the soft X-rays modelling (Table 4). The Fe abundance 
is only slightly super-solar (1.4 ± 0.3 times solar). By includ- 
ing the hard band however, a super-solar Fe abundance is re- 
quired at the 97 per cent level (see results in Table 1 and 4 for 
the hard and soft band respectively). The soft and hard power law 
slopes are consistent with each other with T s — 2.3 ± 0.3 and 
Fh = 2.5 ± 0.3, and the soft power law is required at more 
than the 99.99 per cent level. The hard power law has a normal- 
isation of 5.7^J'g x 10~ 3 ph cm" 2 s" 1 , while the soft one has 
7.0 ± 2.0 x 10~ 5 ph cm -2 s -1 . It should be stressed that the 
soft power law contributes by only ~36 per cent to the 0.5-2 keV 
flux, the remaining ~64 per cent being due to photo-ionized gas 
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Figure 11. The broadband EPIC spectra are shown together with the best- 
fitting spectral model. The model components are as shown in Fig. 6 (bot- 
tom) and Fig. 10 (bottom) for the hard and soft band respectively. 



(~14 per cent) and thermal plasma emission (~50 per cent). The 
possible origins for the soft power law-like continuum compo- 
nent are i) scattered nuclear continuum e.g. from the same photo- 
ionized/collisionally-ionized gas seen in emission in the RGS data, 
or ii) transmitted nuclear continuum (plus reflection) through a par- 
tial covering absorber (see e.g. Levenson et al. 2006 and references 
therein for a discussion). In both cases, we expect less fractional 
variability in the soft than in the hard band because i) the scattering 
medium is most likely extended or ii) the transmitted continuum 
represents 36 per cent only of the soft flux which is dominated by 
extended emission. In both cases the variability would be washed 
out in the soft band much more than in the hard, as observed (see 
Fig 1 and relative discussion). Strictly speaking the first scenario 
would not allow for any short timescale variability in the soft X- 
ray band (and none is detected, see Fig. 1). 

To explore the second solution, we allow the Compton-thin 
absorber to be a partial coverer. The soft power law is removed 
from the spectral model in the attempt of reproducing it with the 
fraction of nuclear emission (hard power law plus reflection) leak- 
ing through the absorber. We find a solution with the same statis- 
tical quality as with the soft power law model (\ 2 = 922 for 820 
dof). The absorber has Nh — 4.1 ± 0.4 x 10 23 cm -2 and covering 
fraction C/ = 98.5 ± 1.1 per cent. However, we also need a fur- 
ther column of 9.5 ± 2.5 x 10 20 cm -2 (in excess of the Galactic 
value) covering totally the nucleus. We conclude that the soft X- 
ray band is dominated by emission lines from both photo-ionized 
gas and thermal plasmas and that the nature of the remaining ^36 
per cent continuum cannot be constrained: a power law continuum 
is statistically equivalent to a ~ 1-2 per cent AGN emission (direct 
continuum plus reflection) leaking through the absorber. 

The observed 0.5-2 keV (2-10 keV) flux is 3.6 ± 1.2 x 
10~ 13 erg cm -2 s _1 (2.1 ± 0.2 x 10~ 12 erg cm" 2 s" 1 ). Our 
best-fitting global model also allows us to compute the absorption- 
corrected luminosity of the different spectral components (see also 
Table 1 and 4). The 2-10 keV observed luminosity is L^ 10 = 
7±2 x 10 42 erg s" 1 only but, if our interpretation of the hard X-ray 
spectrum is correct, the estimated AGN luminosity in the 2-10 keV 
band turns out to be L cst = 4.2 + J:2 x !0 43 erg s" 1 (see Table 1 and 
text). As reported in Table 4, the total 0.5-2 keV luminosity of the 
two thermal components is L th — 1.4 ± 0.4 x 10 41 erg s _1 , while 
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the photo-ionized gas contributes with L ph = 4± 2 x 10 40 erg s _1 
in the same band. Finally, the soft X-ray continuum (scattered or 
leaking through the absorber) has a luminosity of 1.0 ± 0.3 x 
10 41 erg s" 1 in the 0.5-2 keV band. 



8 DISCUSSION 

IRAS 13197-1627 hosts a Compton-thin AGN in which the pri- 
mary X-ray nuclear emission is steep (1^ = 2.5 ± 0.3) and is 
absorbed by a column density of ~ 4 x 10 23 cm -2 of neutral gas. 
A neutral Fe Ka emission line is also clearly detected at 6.4 keV 
and is consistent with being transmitted through the absorber which 
probably has a slightly super-solar Fe abundance. The column den- 
sity inferred from the X-rays is at odds with previous optical stud- 
ies. From the work by Aguero et al (1994) we infer E(B-V)~ 0.47, 
while Cid Fernandes et al (2004) measure Av — 0.38 (note that the 
two values do not satisfy the standard conversion factor Ay-/E(B- 
V)= 3.1). As already pointed out by Maiolino et al (2001a; 2001b) 
E(B-V)/N H ~ 1.2 x 10" 24 , about two orders of magnitude smaller 
than the typical Galactic value (1.7 x 10~ 22 ). This huge discrep- 
ancy is not uncommon in AGNs (Maiolino et al 2001a) and could 
be explained by the presence of large dust grains in the circum- 
nuclear region of AGNs (e.g. due to grain coagulation in dense 
clouds), which results in a featureless extinction curve with reduced 
E(B-V)/iVff . The large grain/coagulation scenario would also sig- 
nificantly reduce the Ay /Nh ratio, as observed here. Other pos- 
sibilities include dust sublimation if the X-ray absorber is suffi- 
ciently close to the luminous nucleus, the idea that most of the 
X-ray absorption is due to the same clouds responsible for the 
broad line region emission, and/or a very low dust-to-gas mass 
ratio for the absorber. Long and short term variability of the X-ray 
absorbing column density would favour explanations based on an 
inner location for the absorber (see e.g. Risaliti et al 2005 and also 
Lamer, Uttley & McHardy 2003). However, previous ASCA and 
BeppoSAX observations show that the flux variability most likely 
occurs at constant spectral slope and absorption. 

The XMM-Newton data also reveal the presence (in both pn 
and MOS detectors) of an absorption line with rest-frame energy 
of 6.81±£ ° l keV, if the redshift of IRAS 13197-1627 is assumed 
(z=0.01654). The absorption line energy points towards an identi- 
fication in terms of an FeXXV Ka line which implies an outflow 
of ~ 5000 km s . However, the outflow velocity is consistent 
within the errors with the systemic receding velocity of the galaxy 
(4959 km s" 1 ) raising the possibility that the gas responsible for 
the absorption is instead local to our own Galaxy or its immediate 
surroundings. 

We find that the most satisfactory description of the hard 
XMM-Newton X-ray spectrum can be obtained by considering 
an additional reflection component from the accretion disc which 
dominates the hard band and also accounts for the 15-30 keV 
Compton hump seen in a previous BeppoSAX observation. The 
reflection-dominated hard X-ray spectrum could be the sign that 
strong gravitational effects are at work in IRAS 13197-1627 (see 
Miniutti & Fabian 2004). Moreover, the disc reflection component 
is correlated with the long-term continuum variability studied by 
using also previous ASCA and BeppoSAX data over 10 years. This 
behaviour is consistent with the light bending model which pre- 
dicts a good correlation between reflection and continuum only 
for reflection-dominated sources, as it appears to be the case here 
but the same prediction holds for any simple reflection model. The 
main advantage of the light bending one is that the large observed 



reflection fraction is automatically explained. It is worth mention- 
ing that, if the extreme observed reflection fraction (R ~ 10) is due 
to light bending effects, the variability behaviour and the value of 
R itself imply that the primary AGN power law continuum is pro- 
duced within 4-5 gravitational radii from the black hole in a com- 
pact and centrally concentrated active region (Miniutti & Fabian 
2004). 

As mentioned by Dadina & Cappi (2004) IRAS 13197-1627 
shares a few properties with the class of Narrow Line Seyfert 1 
(NLS1) galaxies and could represent their obscured counterpart. 
Here we also mention that strong disc reflection components are 
more often seen (or claimed) in NLS 1 than in other type of X-ray 
sources and that the light bending model seems to apply preferen- 
tially to NLS Is, thought to be characterised by high mass-accretion 
rate (e.g. Fabian & Miniutti 2006 for a review). This is not surpris- 
ing given that disc reflection is expected to occur preferentially in 
sources accreting in a radiatively efficient manner, likely to be as- 
sociated with a relatively high accretion rate and with an accretion 
flow extending down close to central black hole where most of the 
accretion power is released. The steep hard X-ray spectra of NLS 1 s 
(and IRAS 13197-1627) may indeed be an indication of a partic- 
ularly high accretion rate in this class (e.g. Shemmer et al 2006). 
As shown in the next Section, a radiatively efficient accretion flow 
seems to be indeed at work in IRAS 13197-1627. In the optical, 
narrow plus broad Balmer lines are detected in IRAS 13197-1627 
with the broad components having a width of ~ 1500 km s" 1 , 
slightly wider than the forbidden lines (~ 700 km s -1 ), suggest- 
ing a Seyfert 1.8 classification but leaving some room for a NLS1 
interpretation (Aguero et al 1994). In this respect it is also very 
interesting to consider that de Robertis et al (1988) reported a num- 
ber of low signal-to-noise broad features in the optical spectrum 
which could be interpreted as Fe II emission, typical of NLS 1 opti- 
cal spectra. 

In the soft band below 2 keV the high-resolution RGS data 
indicate the presence of emission from photo-ionized gas in the 
form of N, O, and Ne emission lines as generally seen in absorbed 
type 2 objects. However, the soft spectrum of IRAS 13197-1627 
is more complex and we detect a plethora of emission lines due 
to a rich Fe L complex from FeXVI and FeXVII, that we interpret 
as due to thermal plasma emission. Both the RGS and pn data can 
be described by considering a two-temperature thermal model in 
which the Fe L emission lines are due to plasmas at ~ 0.6 keV and 
~ 1.4 keV. A power law-like continuum contributes by about 36 
per cent to the soft X-rays and could be due to i) scattering of the 
nuclear emission in the gas/dust revealed in emission/absorption or 
to ii) direct (and reflected) nuclear emission leaking through the 
absorber. The soft X-rays are instead dominated by the thermal 
components (50 per cent), while the photo-ionized gas emission 
represents about 14 per cent of the soft flux. 

8.1 Bolometric luminosity and Eddington ratio 

If the hard X-ray spectrum of IRAS 13197-1627 is indeed dom- 
inated by disc reflection, computing the bolometric luminosity of 
the AGN from the 2-10 keV one is a more difficult and ambigu- 
ous task than usual. The bolometric luminosity of AGNs can obvi- 
ously be computed from the spectral energy distribution and usu- 
ally compares well with the 8-1000/im IR luminosity (Lir). In 
the present case, Lib. ^ 6.7 x 10 44 erg s" 1 (Sanders et al. 2003). 
However about half of it is likely to be associated with a starburst 
(see discussion in next Section) and the AGN IR luminosity can 
be estimated to be L&f N ~ 3.1 x 10 44 erg s _1 . L£? N is most 
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likely produced by dust which intercepts the AGN emission and 
re-radiates it in the IR. In our case, what is the X-ray luminosity 
that will be intercepted by dust and re-radiated in the IR? This is 
not the estimated AGN X-ray luminosity (L|— io> see Table 1 and 
relative discussion) because most of it is bent towards the inner ac- 
cretion disc and never escapes the immediate black hole surround- 
ings. What matters here is the (unabsorbed) observed luminosity 
^2- s io> i- e - me hard X-ray luminosity that can escape the gravita- 
tional sphere of influence of the black hole and makes it to the out- 
ermost nuclear regions irradiating the dust responsible for the IR 
re-emission. Indeed, if one applies standard bolometric corrections 
(Elvis et al 1994) to Lf_\ (see Table 1), one finds a bolometric lu- 
minosity of ~ 2.1 x 10 44 erg s" 1 which is in good agreement with 
Lib., considering the likely large uncertainties in disentangling the 
AGN-starburst contribution to Lir (see next section). 

However, if our spectral modelling and interpretation of the 
hard X-ray spectrum are correct, the above estimate of the bolo- 
metric luminosity has little to do with the real AGN power. Most 
of the power is unobservable directly and is inferred because it pro- 
duces the strong reflection component. This means that L2—10 is 
the quantity that represents the X-ray luminosity extracted from 
the accretion process. Thus, by applying bolometric corrections to 
LI- 10, we can infer the "true" (or better the "estimated") bolo- 
metric luminosity of the AGN, i.e. the luminosity that would be 
observed in absence of strong gravity and light bending effects: 
for consistency, we shall call this quantity Lboi- It turns out that 
Lbol ~ 1-3 x 10 45 erg s" 1 . We stress again that this is not the 
observed bolometric luminosity, but represents instead the lumi- 
nosity that would be observed f the primary X-ray emission was 
isotropic and not gravitationally focused away from the observer. 
Taking our interpretation to the extreme, we then infer that, in ab- 
sence of light bending effects, IRAS 13197-1627 would be seen as 
an ultra-luminous IR galaxy with Lir ~ 10 12 Lq because a much 
higher AGN luminosity would irradiate the dust responsible for the 
AGN IR emission. 

Lfcoj can also be used to compute the Eddington ratio of 
IRAS 13197-1627 if the black hole mass can be estimated. To 
our knowledge there is no black hole mass measurement for 
IRAS 13197-1627. However, the velocity dispersion is cr« = 
143 km s" 1 (Cid Fernandes et al. 2004) and, by using the Mbh-c* 
relation (Gebhardt et al. 2000; Ferrarese & Merritt 2000) as derived 
by Greene & Ho (2006) for local AGN, we estimate a black hole 
mass of 1 -7- 3.5 x 10 7 Mq corresponding to an Eddington luminos- 
ity LEdd = 1.3-^4.6 x 10 45 erg s" 1 . If our interpretation is correct 
and the AGN power is best estimated through L™* h IRAS 13197- 
1627 is radiating between 28 and 100 per cent of its Eddington 
luminosity, although only a fraction of it is actually observed at 
infinity and re-radiated e.g. in the IR. A high Eddington ratio is 
consistent with the detection of a broad Fe line which requires the 
accretion disc to extend down close to the black hole (as likely in 
high Eddington rate sources), with the steep hard X-ray spectrum 
which is a possible indicator of high mass-accretion rate (Shemmer 
et al 2006), and with a NLS 1 association (at least from the X-ray 
point of view). For completeness, we also give the result based on 
the observed 2-10 keV luminosity which is between 4.5 and 16 per 
cent of the Eddington luminosity. All estimates given above should 
be multiplied by about a factor 2 if the highest flux BeppoSAX ob- 
servation is considered. As already mentioned by Dadina & Cappi 
(2004), the BeppoSAX data imply a 2-10 keV AGN estimated lu- 
minosity of ~ 10 44 erg s"\ making IRAS 13197-1627 a type 1.8 
border-line Seyfert/quasar from the X-ray point of view. 



8.2 IRAS 13197-1627: a composite Seyfert/starburst galaxy? 

Besides emission from gas photo-ionized by the nucleus and resid- 
ual/scattered AGN emission, the soft X-ray spectrum is charac- 
terised, and even dominated, by the presence of two thermal com- 
ponents which could be the X-ray signature of star-formation. Cid 
Fernandes et al (2004) estimate that a relatively young (t < 25 Myr 
and possibly t ~ 5 Myr) population of stars represent about the 25 
per cent of the whole population in IRAS 13197-1627. The au- 
thors also point out the detection of the 4680A bump associated 
with Wolf-Rayet stars, an additional indication of star-formation 
activity. 

If the soft X-ray thermal component is indeed due to star- 
formation, its luminosity in the 0.5-2 keV band (1.4 ± 0.4 x 
10 41 erg s" 1 ) can be used to infer the star-formation rate (SFR) 
through the relation SFR X = 2.2 x 10~ 40 Lo.5-2 keV Mq yr -1 
(Ranalli, Comastri & Setti 2003). We do not consider the luminos- 
ity of the photo-ionized gas emission and continuum since these 
are due to AGN emission and we obtain SFRx — 3 1 ± 9 Mq yr - 1 . 
However, some contribution from X-ray binaries may be masked in 
the soft X-ray continuum: we assessed this possibility by forcing 
a r ~ 1.2 power law to be present in the X-ray spectrum (rep- 
resenting the unresolved high mass X-ray binaries contribution, 
see e.g. White et al 1983; Persic et al 2004) obtaining an upper 
limit on any possible additional contribution of 5 Mq yr -1 , i.e. 
SFRx-Slt^Moyr- 1 . 

On the other hand, IRAS 13197-1627 is a luminous infrared 
(LIR) galaxy with Lg-iooo M m = Lir = 6.7 x 10 44 erg s" 1 
(Sanders et al. 2003). The IR luminosity is generally considered 
as a good indicator of the star-formation rate (SFR) in AGN-free 
galaxies, but the reprocessing of the AGN emission in dust can 
dominate the IR luminosity. Rowan-Robinson & Crawford (1989) 
deconvolved the IR spectral energy distribution of IRAS galaxies 
into a disc, Seyfert, and starburst component finding that a ~45 
per cent Seyfert and a ~54 per cent starburst component reproduce 
well the data of IRAS 13197-1627 (with only a marginal disc con- 
tribution). Although the above numbers are probably affected by 
large uncertainties, one can estimate a starburst IR luminosity of 
LfR ~ 3.6 x 10 44 erg s _1 . By assuming this is correct, and by 
using SFRm = 4.5 x lO^Lf^ M Q yr" 1 (Kennicutt 1989), we 
derive SFRir ~ 16 Mq yr -1 , slightly lower than the SFRx esti- 
mate given above. 

As mentioned, IRAS 13197-1627 is also a radio source with 
a Radio luminosity L1.4 ghz ~ 1.6 x 10 30 erg s _1 Hz -1 (Con- 
don et al. 1996). Most of the radio emission is probably due to the 
AGN and, if so, it is not a useful indicator of SFR. At 2.3 GHz, 
the radio core contribution to the total radio flux is about 15 per 
cent (Roy et al 1998). However, even accounting for the radio 
core, IRAS 13197-1627 is still a member of the "radio-excess 
class" among IRAS galaxies (Yun et al 2001), meaning that its 
radio luminosity is in excess from that expected from the well 
known tight radio-IR correlation for normal and starburst galax- 
ies. If we assume that IRAS 13197-1627 would lie on the corre- 
lation if the AGN was removed, the fractional AGN contribution 
in the radio band must be higher than its contribution in the IR. It 
seams reasonable to assume that the AGN contributes by at least 
50 per cent at 1.4 GHz. If so, the AGN-free radio luminosity is 

L ia GHz° C < °- 8 x lo3 ° er 8 s_1 Hz_1 and > fe y aPPtying the re- 
lation SFR1.4 ghz = 0.25 x 10~ 28 L*°^ ree (e.g. Ranalli et al 
2003) we thus obtain SFR1.4 ghz < 20 Mq yr -1 . 

Summarising, we obtain SFRx ~ 31±g 4 Mq yr~\ SFR ffi ~ 
16 Mq yr" 1 , and SFR1.4 ghz < 20 Mq yr" 1 . At face value, 
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the X-ray estimate seems to overestimate the SFR obtained from 
IR and radio. This could be easily explained if a fraction of the 
thermal emission we detect in the soft X-rays was due to gas 
shocked by the radio jet and not to star-forming regions. In the 
prototypical starburst galaxy M82, the plasma is observed to have a 
multi-temperature distribution around 0.4-0.9 keV (Stevens, Read 
& Bravo-Guerrero 2003), similar to the lower-temperature com- 
ponent we detect in IRAS 13197-1627. It is thus possible that the 
high-temperature ~1.4 keV component we detect is not associated 
with the starburst. In this respect, it is interesting to note that if the 
high-temperature X-ray plasma is associated with gas shocked by 
the radio jet and is excluded from the SFRx estimate, we obtain 
SFRx — 18^5°, more in line with the IR and radio estimates. A 
SFR of 15 -r 20 Mq yr _1 probably represents a reasonable value 
which can reconcile the multi-wavelength data. The presence of a 
starburst could be the sign that the nuclear region is gas rich, thus 
providing a large gas supply for the accreting black hole which is 
consistent with the high Eddington ratio inferred above. 
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